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Abstract. Recently, Lobb and Nijhoff initiated the study of variational (Lagrangian) 
structure of discrete integrable systems from the perspective of multi-dimensional con- 
sistency. In the present work, we follow this line of research and develop a Lagrangian 
theory of integrable one-dimensional systems. We give a complete solution of the fol- 
lowing problem: one looks for a function of several variables (interpreted as multi-time) 
which delivers critical points to the action functionals obtained by integrating a La- 
grangian 1-form along any smooth curve in the multi-time. The Lagrangian 1-form is 
supposed to depend on the first jet of the sought-after function. We derive the cor- 
responding multi-time Euler-Lagrange equations and show that, under the multi-time 
Legendre transform, they are equivalent to a system of commuting Hamiltonian flows. 
Involutivity of the Hamilton functions turns out to be equivalent to closeness of the La- 
grangian 1-form on solutions of the multi-time Euler-Lagrange equations. In the discrete 
time context, the analogous extremal property turns out to be characteristic for systems 
of commuting symplectic maps. For one-parameter families of commuting symplectic 
maps (Backlund transformations), we show that their spectrality property, introduced 
by Kuznetsov and Sklyanin, is equivalent to the property of the Lagrangian 1-form to be 
closed on solutions of the multi-time Euler-Lagrange equations, and propose a procedure 
of constructing Lax representations starting from the maps themselves. 

1. Introduction 

It is by now well accepted that considering discrete systems brought a lot of new insight 
and a great deal of simplification into the general theory of integrable systems. This is 
most visible in the field of discrete differential geometry [3J, but also on a more general 
side, a new understanding of integrability of discrete systems as their multi-dimensional 
consistency has been a major breakthrough [2], [To]. This led to classification results [l] 
about discrete 2-dimensional integrable systems (ABS list) which turned out to be rather 
influential. According to the concept of multi-dimensional consistency, integrable two- 
dimensional systems can be imposed in a consistent way on all two-dimensional sublattices 
of a lattices Z m of arbitrary dimension. This means that the resulting multi-dimensional 
system possesses solutions whose restrictions to any two-dimensional sublattice are generic 
solutions of the corresponding two-dimensional system. To put this idea differently, one 
can impose the two-dimensional equations on any quad-surface in Z m (i.e., a surface com- 
posed of elementary squares), and transfer solutions from one such surface to another one, 
if they are related by a sequence of local moves, each one involving one three-dimensional 
cube, like the moves shown of Fig. [TJ 
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Figure 1. Local moves of quad-surfaces involving one three-dimensional cube 



A further fundamental conceptual development was initiated by Lobb and Nijhoff (6) 
and deals with variational (Lagrangian) formulation of discrete integrable systems. So- 
lutions of any ABS equation on any quad surface are critical points of a certain action 
functional obtained by integration of a suitable discrete Lagrangian 2- form (i.e., a skew- 
symmetric real- valued function of an oriented elementary square). Lobb and Nijhoff ob- 
served that the value of the action functional remains invariant under local changes of the 
underlying quad-surface, and suggested to consider this as a defining feature of integra- 
bility. Their results, found on the case-by-case basis for some equations of the ABS list, 
have been extended to the whole list and given a more conceptual proof in [2], and have 
been subsequently generalized in various directions: for multi-field two-dimensional sys- 
tems [7], for dKP, the fundamental three-dimensional discrete integrable system (81, and, 
most recently, for the discrete time Calogero-Moser system, an important one-dimensional 



integrable system 14 



It is the latter generalization, concerned with the Lagrangian structure of integrable 

was 



systems of classical mechanics, that we investigate in the present paper. While 14 



mainly occupied with one concrete example, our aim here is to build of a general theory, 
both in the continuous and discrete time. In (14 , an attempt of such a general theory for 



two-dimensional time was undertaken but unfortunately spoiled by a mistake^ 

In sections [2j |3j we develop the theory for any dimension of the time. The starting 
point is the following problem: under what conditions does a function x of the multi-time 



deliver critical points to functionals 



„ m 

S r = y~]L a (x,x tl , . . . fXt^dtg 



for any smooth curve T in IR m ? At first glance, this seems to be a rather exotic problem. 
However (I am grateful to U. Pinkall and F. Pedit for the corresponding hint), a closely 
related problem lies behind an important theory of pluriharmonic functions and maps. 
Indeed, a pluriharmonic function of several complex variables minimizes the Dirichlet 
functional along any holomorphic curve in its domain. Not surprisingly, both our theory 
and the theory of pluriharmonic functions have one feature in common, namely their rela- 
tion to integrable systems. We derive Euler-Lagrange equations for our problem and show 
that under the Legendre transform they turn into a system of commuting Hamiltonian 
flows with Hamilton functions H a . Moreover, the extremal property formulated above is 
characteristic for commuting Hamiltonian flows (i.e., for an arbitrary system of commut- 
ing Hamiltonian flows one can construct a Lagrangian 1-form with the above mentioned 



1 in their equation (8.116), the derivatives dt%i 'ds and dt^/ds should be interchanged; as a consequence, 
their Euler-Lagrange equation (8.16) is incorrect. 
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extremal property). Note that Hamiltonian flows commute if and only if the Poisson 
brackets of their Hamilton functions are constant, {H a , Hp} = h a p = const. Thus, this 
situation is not necessarily related to complete integrability in the Liouville- Arnold sense, 
where all Poisson brackets should vanish (see, e.g., [TT] for a study of sets of Hamilton 
functions with constant Poisson brackets). It turns out that vanishing of the Poisson 
brackets h a p = is equivalent to the property of the Lagrangian 1-form to be closed on 
solutions of the Euler-Lagrange equations. In other words, for Liouville- Arnold integrable 
systems the value of action Sr on the solutions is invariant under the local changes of the 
curve T. 

In sections |4| [5j we develop the discrete counterpart of this theory. As usual, it turns 
out to be much more transparent and technically simple. Again, the extremal property 
for an arbitrary discrete curve in Z m turns out to be characteristic for a system of m com- 
muting symplectic maps. For Backlund transformations, i.e., for one-parameter families 
of commuting symplectic maps, we show that existence of a certain generating function of 
common integrals of motion, given by the derivative of the Lagrange function with respect 
to the family parameter, is equivalent to closeness of the discrete Lagrangian 1-form. This 
gives a very nice interpretation and explanation of the spectrality property of Backlund 
transformations discovered by Kuznetsov and Sklyanin in [5] . On a more general note, it 
seems to be an important and difficult problem to find general enough conditions assuring 
integrability of a given system of commuting symplectic maps. 

For Backlund transformations, we push forward the idea that such a standard integra- 
bility attribute as Lax or zero curvature representation (yielding, among other features, 
existence of a sufficient number of integrals of motion) can be extracted from the maps 
themselves. To put this more aphoristically, a family of Backlund transformations serves 
as its own Lax representation. This idea is conformal with the corresponding idea for dis- 
crete two-dimensional systems with the property of multi-dimensional consistency, which 
also serve as zero curvature representations for themselves (cf. [2j, [lO]). In section |6j we 
demonstrate, by the way of example, how our theory works for Backlund transformations 
for the classical Toda lattice. 

2. Basic construction in the continuous time case 

Consider a 1-form on IR m (which we call the mult i- dimensional time or simply multi- 
time), whose coefficients are functions on X x X m : 

m 

C(x,vi, ...,v m ) = y^Lajx^i, . . .,v m )dt a . (1) 

a=l 

Here X is the configuration space of a variational problem (usually we take X = M. N ). 
We are looking for a function x : lR m — > X which solves the following variational problem. 
For any (smooth) curve T : [0, 1] — » lR m , it is required that x delivers a critical point for 
the functional 

= J C(x,x tl ,...,x tm ). (2) 
We will prove the following results. 

Theorem 1. The function x : M. m — > X delivers a critical point for the functional Sr 
for any smooth curve T, if and only if the following conditions (which we call multi-time 
Euler-Lagrange equations,) are satisfied: 



YURI B. sums 

For any 1 < a 7^ < m, 



^^(x,x tl ,...,x tm ) = 0; (3) 

• The following relations hold true: 

dL 1 dL m 

— (x,x tl ,...,x t J = ... = - — (x,x tl ,...,x tm . 4 

VKe denote the common value of these functions by p : R m — > X ; 

• This function p satisfies the following differential equations: 

dp dL a 

^- = ^(x,x tl ,...,x t J. (5) 

The system of multi-time Euler-Lagrange equations is usually overdetermined, and it 
is not easy to analyze its compatibility in general. Indeed: 

• Among m{m— 1) + (m — 1) = m 2 — 1 (vector) equations (|3]), Q, there should be 
precisely m — 1 independent ones. 

• Differential equations ^ are second order partial differential equations, and the 
same is true for their compatibility conditions, 

d ( dLp \ d ( dL r 



dta v dx {x, x tl , . . . , x t J) - — {—(x, x tl , . . . , x tm )j = 0. (6) 

One can expect that these compatibility conditions should be satisfied by virtue 
of equations ^ themselves. 

Example. Consider the two-time Lagrangian 1-form with the following components: 

N N 



Li(x,x tl ) = I^((x fe ) tl ) 2 -X>^-*<% (7) 

k=l k=l 
N ^ N N 

L 2 {x,x tl ,x t2 ) = ^(a;fc) tl (x fc )t 2 --^((a; fe ) tl ) 3 -^e a;fc + 1 - Xfc ((x fc )i 1 + (a; fc+ i) tl ). (8) 
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k=l k=l k=l 



Then one of the two equations (|3j) is trivially satisfied, while the second one reads. 

{x k )t 2 = {{xk)nf + e x ^~ Xh + e x ^~\ (9) 

Equation |4]) is satisfied identically. Finally, equations ^ read: 

(x k )t ltl = e x ^~ x * - e x *-**~\ (10) 
(x k )t lt2 = e a;fc + 1 - a;fc ((a; fc ) tl + (x fc+1 ) tl )-e :Efc - a;fe - 1 ((x fc „ 1 ) tl + (x fc ) tl ). (11) 

Equations (J9]) ( 1 1 ) constitute the two-time Euler-Lagrange equations for the above La- 
grangian 1-form. One shows by a direct computation that the compatibility condition of 
the latter two equations, ((xk)t 1 t 1 )t 2 — ((x k ) tlt2 ) tl , is satisfied by virtue of equations ^ 
and (10). A well-posed initial value problem for equations (J9]) ( 1 1 ) can be obtained by 
prescribing x and x tl at one point (ti,t 2 ) = (0,0) of the two-time. 

In order to arrive at more tractable conditions, we adopt the following definition. 
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Definition 1. We say that the multi-time Lagrangian 1-form £ is Legendre transformable, 
if the set of equations consisting of Q) and 

9L± . . dL m 

— {x 1 x tll ... 1 x tm ) = ... = -^{x,x tl ,...,x tm ) = p (12) 

can be solved for x ta in terms ofx,p: 

x t a = x ta (x,p), ct = l,...,m. (13) 

Clearly, for Legendre transformable 1-forms the right-hand sides of differential equa- 
tions ([5]) can be expressed through X ? J)j clS well. 

Example, continued. In the example above, equations (13) take the form 

(ac*)ti = Pk, (14) 
{x k )t 2 = pl + e x ^- x " +e Xk ~ Xk -\ (15) 
while equations ^ can be put as 

(Pk) tl = e x ^~ Xk - e"-\ (16) 
(Pk)t 2 = e x ^- Xk ( Pk +p k+1 )-e Xk - x ^(p k ^+ Pk ). (17) 
One easily sees that equations ( 14 ) , ( fi~6| are Hamiltonian with the Hamilton function 

N N 
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k=l k=l 



while equations (15), (17) are Hamiltonian with the Hamilton function 

^ N N 

k=l k=l 

These observations illustrate the following general statement. 

Theorem 2. Suppose that the multi-time Lagrangian 1-form C is Legendre transformable. 
Define the Hamilton functions 

H a (x,p) = {p,x ta )-L a (x,xt 1 ,...,x tm ), a = l,...,m (18) 

(it is understood that all x t „ on the right-hand side of the latter formula are expressed 
according to (13)). The functions x,p : M m — > X satisfy Hamiltonian equations of motion: 

dx = dH a dp dH a 

dt a dp ' dt a dx 

Hamiltonian flows with Hamilton functions H a commute, so that their pairwise Poisson 
brackets are constant: 

{H a , Hp} = h a p = const. (20) 

Theorem 3. The following identities hold true on solutions of the Euler-Lagrange equa- 
tions: 

d d 

Qj-Lp(x, x tl ,..., x tm ) — — L a (x, x tl ,..., x tm ) = h a p, (21) 

with the constants h a p from (20). In particular, if the Hamilton functions H a are in 
involution, h a p = 0, then the form C is closed on solutions of the Euler-Lagrange equa- 
tions, so that the action functional Sp in M) does not depend on the choice of the curve 
T connecting two given points in W m . 
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Conversely, any system of commuting Hamiltonian flows admits a Lagrangian formu- 
lation with a 1-form satisfying conditions of Theorem [T] A trivial version of such a 
formulation is obtained upon individual Legendre transformations of all functions H a , so 
that each L a only depends on its own velocity: L a = L a (x,x ta ). However, this construc- 
tion might be either impossible (if some of Legendre transformations are non- invert ible, 
i.e., relations Xt a = dH a /dp are not solvable for p in terms of x,Xt a ), or just imprac- 
tical (if the above mentioned relations are solvable for p but the resulting expressions 
are inconvenient, i.e., given by algebraic functions of high degree). The following theorem 
provides us with a convenient construction which works whenever just one of the Legendre 
transformations is easily invertible, which is usually the case in applications. 

Theorem 4. For any system ofm commuting Hamiltonian flows with Hamilton functions 
H a (x,p), a = 1, ... ,m, one can find a multi-time Lagrangian 1-form satisfying conditions 
(and statements) of Theorem^ In particular, if the Hesse matrix d 2 Hi/dp 2 is non- 
degenerate, so that the relations x tl = dH\/dp can be solved for p = p(x,x tl ), then one 
can chose the components of £ as follows: 

Li(x,x tl ) = (p,x tl ) - Hi(x,p), (22) 
L a (x,x tl ,x ta ) = (p, x ta ) - H a (x,p), a = 2,...,m (23) 

(with the understanding that p on the right-hand sides of all these formulas is expressed 
through x,x tl ). 

3. Proof of main theorems in the continuous time case 

Proof of Theorem^ For a function x : R m — > X and for a fixed curve T : [0, 1] — > M. m , 
we denote by xr = x o V : [0, 1] — > X the corresponding curve in X. A variation of the 
action functional S-p due to an arbitrary variation 5x : IR m — > X of the function x can be 
written as 

** = rt((^- fa ) + g<fe-^)) I '-W' fa < 24 > 

Here (and in similar context below) it is understood that all terms in the parentheses are 
evaluated at T(s). The partial derivatives of Sx, when restricted to T, satisfy the following 
two conditions: 

• The derivative of 5x in the tangential direction of T is given by 

m 

^2T t a (8)6x ta = (6x r )' i (25) 

a=l 

where d~xr '■ [0, 1] — > X is an arbitrary variation of the function xr- 

• The derivatives of Sx in the normal directions to the curve T are arbitrary, i.e., 
are independent of Sxp. More precisely, if nx(s), . . . , n m _i(s) is an arbitrary basis 
of (r'(s))- 1 -, the normal space to the curve T at the point T(s), then 

m 

^riypSxtp = 5y 71 7 = 1,..., 771 - 1, (26) 

where n 7/ g are the components of the vector n 7 , and 8y 1 are m — 1 arbitrary 
functions on V. 
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The solution of system ( 25 ) , ( 26 ) is most easily found, if we assume vectors n\ (s) , . . . , n m _i (s) 
to be orthonormal, and is then given by 

r (?) m_1 

Sxta = W(sW {5xr) ' + Ji n ~< a{s)5y "- (27) 



Now we are in a position to find necessary and sufficient conditions for x to ensure 5 St = 0. 
We have: 



5St 



i / m 



E^ 



s)^tt— , 5xr ) ds 



dx 

m 

+E/ <EE^: 



a=l 

m—l r \ I m m 



a=l 0=1 



|F(s)||2 



7=1 



[s)^-^-n ry p{s),5y 1 ) ds. 



o 



<9t> 



(28) 

or : / 

K a=l /3=1 

By the usual integration by part argument in conjunction with the fact that Sxt and 5y^ 
(7 = 1, . . . ,m — 1) are arbitrary independent functions on T, we come to the following 
necessary and sufficient conditions of 5 St = 0: 

ri,Wr» di c 



E r « 

a=l 



0L a d ^ — > ^ — > 

dx ds ' ' 



a=l 0=1 



|P( S )||2 fyg 



0. 



<9L 



(29) 



(30) 



a=l /3=1 

The arguments of the functions dL a /dx, dL a /dvp in these formulas are (x,x tl , . . . ,x tm ) 
for x : lR m — > X, evaluated at T(s). Once again, equations (29), (30) characterize functions 
x delivering a critical point for the functional St for a fixed curve T. 



Since the curve T is arbitrary, we can further argue as follows. Equation (30) says that 

m m r, T 



dv 







a=l (3=1 

for an arbitrary pair of orthogonal vectors T'(s) and n(s) at an arbitrary point T(s) G M. m . 
A necessary and sufficient condition for this is given by (J3|, Q. Substituting the latter 
into equation ( 29 ) results in 



En, 

a=l 



^ dL a dp 
dx ds 



0. 



which is equivalent to (J5|. 



Proof of Theorem This is a standard computation, not much different from the one 
in the 1-time situation. We compute: 

' dL a dx t 



dH a 
dp 



dx t 
dp 



y 

^\dvp' dp 



Indeed, due to (J3|, the sum reduces to one term with (3 = a, which cancels with the 
second summand on the right-hand side, due to Q. This proves the first equation of 
motion in (19). 
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Literally the same argument justifies the following computation: 

dH a I dx ta \ dL a ^ I dLa dx tfS \ dL a 

0=1 



dx \ ' dx I dx ^— ' \ dvp 1 dx / dx 



This proves the second equation of motion in (19), by virtue of ([5]): 

dp dL a dH a 
dt a dx dx 

Thus, Hamiltonian equations of motion are verified. Commutativity of flows is a mere 



reformulation of the existence of functions x,p : lR m — > X which solve all equations (19) 
simultaneously. ■ 

Proof of Theorem [3| A very similar computation: 

9 t ( rr r )-/ dL " dx\ ^/dL a dx h 

Upon using equations of motion (JiJ), we come to 

d , / dp dx \ / d 2 x 

-L a {x, x tl ,..., x tm ) = < — — , — — ) + < p, 



dtp " v ' l1 ' ' tm/ \dt a 'dtf}/ \'dt a dtf) 
This yields 

d , . d . . / dp dx\ I dp dx 

—L a {x,x tl ,. ..,x t J- — L p [x, x ^,...,x u ) = ^—,—j-^—,— 



and, by virtue of Hamiltonian equations of motion (19) 



d _ , , d ldE a dHA /dHp dH Q 

—L a (x, x tl , . . .,x tm ) - Wa Lp(x,x tl , . . . , x tm ) = - ^ — , + — 

which is the desired result. ■ 

Proof of Theorem [^} A direct verification of conditions of Theorem [TJ ■ 



4. Basic construction in the discrete time case 

In this section, we consider a sort of a discretization of the previous construction. We 
are looking for functions x : Z m — > X — M. N which deliver critical points for the following 
variational problem. Let F be an arbitrary discrete curve (path) which is a concatenation 
of a sequence of directed edges in Z m such that the endpoint of any t k is the beginning 
of efc+i, we set 

S T = Y J ^k)- (31) 

fcez 

Here £ is a discrete 1-form, i.e., a skew-symmetric function on directed edges of the regular 
square lattice Z m , defined as follows: 

C(n, n + e a ) = A a (x, x a ) C(n, n — e a ) = — A a (x_ a , x). 
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Here A a : X x X — > K are local Lagrangian functions corresponding to the edges of 
the a-th coordinate direction, and the following abbreviations are used: x for x(n) at a 
generic point n G Z m , and then 

x a = x(n + e a ), x_ a = x(n — e a ), a = l,...,m. (32) 

Here e a stands for the unit vector of the a-th coordinate direction. 

Theorem 5. The function x : Z m — > X delivers a critical point for the functional S-p 
for any discrete curve T, if and only if the following conditions (which we call multi-time 
discrete Euler-Lagrange equations j are satisfied: 

dA a (x_ a ,x) | dApfaxp) = Q 3 



dx dx 
dA a (x- a ,x) dAp(x-p,x 



dx dx 
dA a (x,x a ) dAp(x,xp) 



0, (34) 
0. (35) 



dx dx 

In other words, there exists a function p : Z m — > X satisfying all the relations 

dA a (x, Xa) tnr\ 

P = , a = l,...,m, (36) 

9A a [X-. a ,X) /nn\ 

V = q- , a = l,...,m. (37) 

Continuing analogy with the continuous time case, we introduce the following defini- 
tion. 

Definition 2. We say that the multi-time discrete Lagrangian 1-form C is Legendre 
transformable, if all the equations (36) can be solved for x a in terms ofx,p. 

Theorem 6. Suppose that the multi-time discrete Lagrangian 1-form is Legendre trans- 
formable. Then equations 

dA a (x,x a j dA a [x^x Q ) /oo\ 

V = 7, , Pa = « , (38) 

dx dx a 

determine a symplectic map F a : (x,p) (->■ (x a ,p a ). These maps for different a commute: 

F a o F p = Fp o F a . (39) 

Theorem 7. The following identities hold true on solutions of the discrete Euler-Lagrange 
equations: 

A a (x, x a ) + Ap{x m x a p) - A Q (x (3 , x a p) - Ap(x, xp) = i a p = const. (40) 

In particular, if all these constants E a p vanish, then the discrete 1-form C is closed on 
solutions of the Euler-Lagrange equations, so that the action functional Sr does not depend 
on the choice of the curve F connecting two given points in Z m . 

In the full generality, the above discrete theory does not lead to any statements about 
integrability of the maps F a . However, such statements become possible for an important 
class of examples, namely, for Backlund transformations. We understand Backlund trans- 
formations as a one-parameter family of commuting symplectic maps. Thus, our point of 
view is in a sense opposite to that of Kuznetsov and Sklyanin in [5J. While the primary 
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feature of Backlund transformations for them was the existence of a common complete 
set of integrals in involution coming from a common Lax matrix (and commutativity was 
considered as a consequence of this property by virtue of the discrete Liouville- Arnold 
theorem), we propose to put an emphasis on the commutativity property. Actually, we 
even do not pre-suppose the existence of the Lax representation for the Backlund trans- 
formations, but rather derive it from the maps themselves. 

Thus, we consider the family of maps Fx : M. 2N — > M? N , (x,p) i— > {x,p), given by 



equations of the type (38): 



dA(x,x;X) _ dA(x,x;X) 

F\: p= q , p = 7^ • 41 

ox ox 

When considering a second such map, say F^, we will denote its action by a hat: 

dA{x,x;p) _ dA{x,x;p) 

' P= dx ' P= dx— (42) 

Assuming that F\, F^ commute for any A and p, we write the result of Theorem [7] as 

A(x, x; A) + A(x, x; p) — A(x, x; p) — A(x, x; A) = £(X, p). (43) 

Clearly, the possible dependence of this constant on the parameters A, p is skew-symmetric: 
£{X,p) = -£{p,X). 

Theorem 8. For a family of Backlund transformations, the discrete 1-form C is closed 
on solutions of the multi-time Euler Lagrange equations, i.e., £(X, p) = 0, if and only if 
dA(x, x; X)/dX is a common integral of motion for all F^. 

The latter property is a re-formulation of the spectrality property of Backlund transfor- 
mations introduced by Kuznetsov and Sklyanin [5] but stripped of its mystical flavor by 
not mentioning the additional structure of the zero-curvature, or Lax, representation. 

Turning the subject of the Lax representation, we would like to pursue the one- 
dimensional analogue of the basic idea of consistency as integrability j2], [10]. According 
to this idea, a discrete multi-dimensionally consistent two-dimensional system serves as 
its own zero curvature representation. Our message concerning the one-dimensional case 
is: 

Given a one-parameter family of commuting symplectic maps, one can derive a 
Lax (or zero- curvature) representation with a spectral parameter from the maps 
themselves. 

This can be considered as a further step to de-mystifying integrable systems. Indeed, 
commutativity is a simple property allowing for a direct (may be, computer aided) check, 
while Lax representations are transcendental objects whose mere existence has always 
been considered as the most mysterious feature of integrable systems. 

We will not give a general derivation, but instead present it in the particular case of 
Backlund transformations for the Toda lattice described in Section [6| An extensive body 
of concrete results for further families of Backlund transformations will be presented 
elsewhere. 

5. Proofs of main theorems in the discrete time case 

Proof of Theorem [5| At any point of a discrete curve there meet two directed edges. 
It is easy to understand that, up to an overall change of orientation of the curve and up 
to permutations of indices, only four possibilities exist, depicted on Figure [2} 
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X/3 



X-r 



X 



X 



X f 



X. 



X 



(b) 



X 



x, 



(d) 



X-0 



Figure 2. Four type of vertices of a discrete curve. Cases (a) and (b): two 
equally (positively or negatively) directed edges meet at n, of one and of two 
different coordinate directions, repectively. Case (c): a positively directed 
adges followed by a negatively directed edge. Case (d): a negatively directed 
edge followed by a positively directed edge. 



Correspondingly, there are four types of the Euler-Lagrange equations. The first one 
is for a straight piece of a discrete curve, as on Fig. |2j (a): 

OA ) OA 



+ 



dx dx ' 

while the other three correspond to the corner type pieces of a discrete curve, as on Fig. |2j 
(b)-(d), and are given in (33)-(35). It is important to observe that these equations are 
not independent, for instance (44) follows from (33) and (35). Existence of the function 
p : Z m — > X is a direct corollary of equations (33)-(35). ■ 



Proof of Theorem [6[ The symplectic property of the maps F a was proven by Moser- 
Veselov M. The commutativity of maps F a is a mere restatement of the existence of 
functions x,p : Z m — > X which solve all equations (36), 
that a function (x,p) : Z m — y M? N such that (x(n + e c 
defined by a single initial value (x(0),p(0)). ■ 



37j) simultaneously. We remark 
,p(n + e a )) = F a {x{n),p{n)) is 



Proof of Theorem^ Denote the left-hand side of equation (40) by Vt. It is a function 
on the manifold of solutions of discrete Euler-Lagrange equations, which is, according to 
the last remark in the previous proof, (2iV)-dimensional, as it can parametrized by (x,p), 
or by (x, x a ) with some fixed a. It is enough to prove that dQ/dx = and dQ/dx a = 0. 
We prove a stronger statement: if one considers Q as a function on the (4iV)-dimensional 
manifold, parametrized by x, x a , xp, and x a p, then the gradient of this function vanishes 
on the (2iV)-dimensional submanifold of solutions of discrete Euler-Lagrange equations. 
But this is obvious, since vanishing of the partial derivatives of Q with respect to its 
4 vector-valued arguments is nothing but the discrete Euler-Lagrange equations at the 
corresponding points. ■ 



Proof of Theorem^ Due to skew-symmetry, £(X, fi) = is equivalent to d£(X, p)/d\ 
0, that is, to 

<9A(x,x;A) dA(x, x; X) 



dX 



dX 







see equation (43)). This is equivalent to saying that dA(x, x; X)/dX is an integral of 



motion for F M . 
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6. Example: Backlund transformations for Toda lattice 

Here we illustrate the main constructions of the discrete time case by the well-known 
example of Backlund transformations for the Toda lattice 13 , [5] , also known as discrete 
time Toda lattice [12]. From this one-parameter family of commuting symplectic maps, 
we consider just two (for notational simplicity only; considering any m would go along 
literally the same lines). Our maps F\,F^ : M. 2N — > M. 2N are given by equations of the 
typeg: 

F x : p k = - [e Xk ~ Xk - 1) + Xe*"-*"- 1 , p k = - (e Xk ~ Xk - l) + Xe Xk+1 ~ Xk , (45) 
A A 

Fn : p k = - (e Xk ~ Xk - 1) + fie Xk ~ Xk -\ p k = - ( e Xk ~ Xk - l) + ^e Xk +^~ Xk . (46) 

The corresponding Lagrangians are given by 

N N 



A{x,x;X) = jJ2(e Xk - Xk -l-(x k ~x k ))-Xj2e Xk ^- Xk , (47) 

k=l k=l 

N N 

A(x,x;n) = -^{e Xk ~ Xk -\-{x k -x k ))-^ eXk+1 ~" k - ( 48 ) 



r k=l k =l 

These maps are considered under open-end boundary conditions (x = oo, xn + i = — oo) 
or under periodic boundary conditions (all indices are taken modiV, so that x = x^, 
Xn+i — X\). The corresponding discrete Lagrangian 1-form is Legendre transformable. 
To see this, consider the map F\. In the open-end case, the first equations in (45) are 
uniquely solved for x%, x 2 , ■ ■ ., xn (in this order), to give 

\2pX2~Xl 

e Zl-*l = l+\p u e *2-X2 =1+Xp7 



1 + Api 

^2 e x N -x N _ 1 

y2pX N _ 1 -x N _ 2 

1 + XpN-2 — . \^ e X2-xi 



= 1 + A ^ >2^-x„_, • (49) 

1 + Xp N - 



1+Xpi 

In the periodic case e Xk ~ Xk can be expressed as analogous infinite periodic continued 
fractions and are, therefore, double-valued functions of (x,p). 

It is well known that the maps Fa, F^ commute, see, e.g., [13], [5], [12] (however, the 
very notion of commutativity for double-valued maps, relevant for the periodic case, re- 
quires for a more careful discussion, which we will provide in a separate publication). The 
most direct way towards the proof of commutativity relies on the following superposition 
formula: 

» ~ ~ \ p Xk+l _ n p x k+l 

pXk-Xk-Xk+Xk+i = {^Q\ 
\ e x k _ ^Xk ' 

which is also equivalent to either of the two formulas 

^ e X k + l~Xk + l _ J^gZfc + l-Zfc + l _ \pX h -Xk _|_ \pXk-X k _ Q (51) 

A \x A n 

X e ^k + i-xk _ ^k+i-xk _ \ e x k +i-z k + f^h+i-xh = o. (52) 
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Thus, for any initial point (x(0),p(0)) maps F\ and F^ determine a function x : Z 2 — > M. N 
which, according to Theorem [5j delivers a critical point to the action defined by an 
arbitrary curve T in Z 2 and by the Lagrangian 1-form C with the Lagrangians A(x, x; A) 
and A(x,x;/i) on the edges parallel to the first (resp. second) coordinate direction and 
satisfies the following discrete Euler-Lagrange equations: 



Pk 



^ x k— x k l) + Xe Xh ~ a; * -1 



(e* k ~ Xk - 1) + fie 



*^fc *^fc — 1 



(< 



x k -x k 



1) + Ae^ +1 " 



x A . 



(e 8 *-** - 1) + /ie 



(53) 
(54) 



/x /i 
Moreover, we can prove that the Lagrangian 1-form C is closed on solutions. For this, we 



first show that this property, i.e., £(X, fi) = for £(X,fi) from (43), is equivalent to 



N 

D 

k=i 



x k - x k -x k + x k ) = 



N 

n 

k=l 



,x k -x k -x k +x k 



1. 



(55) 



For this aim, we observe that, by virtue of the superposition formulas (51), (52), most of 



the terms on left-hand side of (43) cancel, leaving us with 

(\ l\ N ^ 
£(X,fi) = I y^(%k ~x k -x k + x k ). 



Alternatively, one can refer to Theorem |8j Indeed, one easily computes: 

dA(x, x; A) 1 



A I N 

k=l k=l 



Xk 



the first sum on the right-hand side being an obvious integral of motion. Now the desired 



result (55) follows in the periodic case by multiplying equations (50) for 1 < k < N, in 



the open-end case equation (50) holds true for 1 < k < N — 1 and has to be supplemented 



by the boundary counterparts 



Ae 



fie 



Xjy" — Xjy — Xjy 



A 



(56) 



A — fi Xe XN — fie XN 

Thus, spectrality of Backlund transformations for the Toda lattice is a direct consequence 



of the superposition formula (50) 



Finally, we turn to the topic of Lax, or zero-curvature representations for Backlund 
transformations. We consider the map F\ as the basic object (it represents actually 
the whole family of maps, since the parameter A is arbitrary), and will produce the 
zero curvature representation from the map F^, with /i playing the role of the spectral 
parameter. For this aim, we will extract the "wave functions" tpk from the variables x k - 
Indeed, setting 



?x k -x k 



(which defines ipk up to a common factor), we re- write the first equation in (46) as 

1 

M 



+ Pk I tpk - e^-^-^k-i- 
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This is nothing but the first (spectral) part of the inverse spectral formulation of the Toda 
hierarchy. 

This first part of derivation was known from the early days of the soliton theory For 



instance, it can be found in the paper 13 by Wadati and Toda from 1975. They treated 
Backlund transformations as symplectic maps commuting with the respective continuous 
time flows, and correspondingly demonstrated how to derive a Lax representation for 
a continuous flow from a Backlund transformation. This amounts to deriving a linear 
differential equation describing the temporal evolution of the wave function ip k . Thus, 
we depart from their derivation only now, by considering a commuting map instead of a 
commuting flow, which will lead to a difference equation for the discrete time evolution 
of the wave function ip k . It is amazing that this seemingly simple step (leading to a great 
conceptual simplification) took about 40 years to be done. 

Now, we turn to deriving a consequence of the commutativity of F\ and F„, expressed 



via the superposition formula ( 50 ) . It can be re- written as 



H^^- = e k k = = /! — : , 

ip k A — jj / e Xk ~ Xk e Xk ~ Xk Aipk+i — e Xk ~ Xk ipk 

so that the second (evolutionary) part of the inverse spectral formulation is obtained: 

i>k = -A^fc+i + e* k ~ Xk ilj k . 

We get a standard zero-curvature representation of the Backlund transformation for 
the Toda lattice (see, e.g., [H], [12]) upon introducing an auxiliary wave function <f) k = 
—e~ Xk - 1 ipk-i- Then we arrive at 



ipk+i\ (l/v+Pk e Xk \ (ip k 




Ae-^-i 1 



(57) 



(58) 



Denoting the 2x2 matrices appearing in (57), (58), by Lk(fJ>), Vfc(/i), respectively, one 



arrives at the zero curvature representation of the map Fx with the spectral parameter \i: 

L k (ii)V k (ii) = V k+1 {y)L k {n), 
which, in particular, yields that the trace of the monodromy matrix 

T(n) =L N (n) . . . L 2 (fi)L 1 (n) 

is an integral of motion (or, better, a generating function of integrals of motion) of the 
map F\. By the way, in the open-end case one easily finds with the help of continued 
fractions for e Xk ~ Xk given in M9J) that 

N 

Y[ e Xk Xk = A Ar (T(A))n. 

k=X 

In the periodic case, one has to replace the right-hand side by one of the eigenvalues 
of the matrix X N T(X). In both cases, the zero-curvature representation yields that the 
left-hand side is an integral of motion. Recall that this statement is equivalent to the 
discrete Lagrangian 1-form being closed on solutions. 
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7. Conclusions 

In this paper, we developed a variational characterization of systems of commuting 
Hamiltonian flows (resp. of commuting symplectic maps). This characterization can be 
formulated as extremizing the action obtained by integration of a Lagrangian 1-form along 
an arbitrary curve in the multi-dimensional time. There remains a challenge to develop a 
similarly complete theory for two- and three-dimensional systems, both in the continuous 
and discrete contexts. This is the subject of the ongoing research. 

This research is supported by the DFG Collaborative Research Center TRR 109 "Dis- 
cretization in Geometry and Dynamics". 
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